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In their ideal structure, double perovskites like S^FeMoOe have alternating Fe and Mo along 
each cubic axes, and a homogeneous ferromagnetic metallic ground state. Imperfect annealing leads 
to the formation of structural domains. The moments on mislocated Fe atoms that adjoin each 
other across the domain boundary have an antiferromagnetic coupling between them. This leads 
to a peculiar magnetic state, with ferromagnetic domains coupled antiferromagnetically. At short 
distance the system exhibits ferromagnetic correlation while at large lengthscales the net moment is 
strongly suppressed due to inter-domain cancellation. We provide a detailed description of the spin 
wave excitations of this complex magnetic state, obtained within a 1/S expansion, for progressively 
higher degree of mislocation, i.e., antisite disorder. At a given wavevector the magnons propagate at 
multiple energies, related, crudely, to 'domain confined' modes with which they have large overlap. 
We provide a qualitative understanding of the trend observed with growing antisite disorder, and 
contrast these results to the much broader spectrum that one obtains for uncorrelated antisites. 



I. INTRODUCTION 

Double perovskite (DP) materials with general formula 
A2BB'06 have generated a great deal of interest^ both 
in terms of their basic physics as well as the possibility 
of technological applications. In particular, S^FeMoOe 
(SFMO) shows high ferromagnetic T c , ~ 420K, large elec- 
tron spin polarisation (half-metallicity) and significant 
low field magnet or esistance^. 

The ferromagnetic coupling between the S = 5/2 lo- 
calized magnetic moments in SFMO (Fe 3+ ion, 3d 5 state) 
is driven by a "double exchange" mechanism, where elec- 
trons from Mo delocalise over the Mo-O-Fe network. The 
B (Fe) ions order ferromagnetically while the conduction 
electrons that mediate the exchange are aligned opposite 
to the Fe moments, leading to a saturation magnetization 
of AfiB per formula unit in ordered SFMO. However, the 
large entropy gain from disordering promote 'antisite dis- 
order' (ASD) whereby some B ions occupy the positions 
of B' ions and vice versa. 

There is clear evidence now that B-B' mislocations are 
not random but spatially correlated 4 5 . While ASD sup- 
presses long range structural order, electron microscopy 4 ^ 
and X AFUP reveal that a high degree of short range order 
survives. The structural disorder has a direct magnetic 
impact. If two Fe ions adjoin each other the filled shell 
d 5 configuration leads to antiferromagnetic (AFM) su- 
perexchange between them. The result is a pattern of 
structural domains, with each domain internally ferro- 
magnetic (FM) while adjoining domains are AFM with 
respect to each other. This naturally leads to a suppres- 
sion of the bulk magnetisation with growing ASD. 

Domain structure has been inferred in the low doping 
manganites as well, due to competing FM and AFM in- 
teractions. Inelastic neutron scattering in those materials 
suggest the presence of FM domains in a predominantly 
AFM matrix, and allows a rough estimate of the domain 
size^El ^ e aim to provide a similar framework for in- 
terpreting the magnetic state and domain structure in 



the DP from spin wave data. Our main results are the 
following. 

(i) We compute the dynamical magnetic structure fac- 
tor, that encodes magnon energy and damping, within a 
1/S expansion of an effective Heisenberg model chosen 
to fit the electronic model results, (ii) The magnon data 
is reminiscent of the clean limit even at maximum ASD 
(50%), where the bulk magnetisation vanishes due to in- 
terdomain cancellation, (iii) We suggest a rough method 
for inferring the domain size from the magnon data and 
check its consistency with the ASD configurations used, 
(iv) We demonstrate that uncorrelated ASD leads to a 
much greater scattering of magnons and a much broader 
lineshape. This suggests that in addition to XAFS and 
microscopy, neutron scattering would be a sensitive probe 
of the nature of disorder in these materials. 

The paper is organized as follows: In Sec. II we discuss 
the generation of the structural motif, the solution of the 
electronic problem, and the estimation of exchanges for 
an effective Heisenberg model. In Sec. Ill we recapitu- 
late the spin-wave formulation for non collinear phases 
and present the magnon spectrum obtained for the dif- 
ferent disordered configurations. In Sec. IV we discuss 
the results, attempting to analyse the magnon spectrum 
for correlated antisites in terms of confined spin wave 
modes, and contrasting the result to magnons in an 'un- 
correlated' antisite background. 

II. EFFECTIVE MAGNETIC MODEL 

A. Structural motif 

Given the similar location of the B and B' ions (at 
the center of the octahedra) the tendency towards defect 
formation is more pronounced in the DP's. This ten- 
dency of mislocation interplays with the inherent B-B' 
ordering tende ncy and creates a spatially correlated pat- 
tern of antisites^El rather than random mislocation. To 



model this situation we have used a simple "lattice-gas" 
modeP. On proper annealing it will go to a long range 
ordered B, B', B, B'... pattern. We frustrate this by 
using a short annealing time to mimic the situation in 
the real materials. We encode the atomic positions by 
defining a binary variable 77^ such that rji = 1 when a 
site has a B ion, and rji = when a site has a B' ion. 
Thus for an ordered case we will get 77's as 1, 0, 1, 0, 1, 0... 
along each cubic axes. The B-B' patterns that emerge 
on short annealing are characterised by the structural 
order parameter S = 1 — 2x, where x is the fraction of 
B (or B') atoms that are on the wrong sublattice. We 
have choosen four disordered families with increasing dis- 
order for our study. One structural motif each for these 
families is shown in the first column of Fig.l, with pro- 
gressively increasing disorder (from top to bottom). We 
plot g(ri) = (rji — i)e 27r( ^ + ^) as an indicator of struc- 
tural order. For a perfectly ordered structure g(ri) is con- 
stant. The pattern along the first column are different 
realisations of ASD with S = 0.98, 0.88, 0.59, 0.17 (top to 
bottom). We solve the electronic-magnetic problem on 
these structural motifs. 



B. Electronic Hamiltonian 

To study the magnetic order we use the electronic- 
magnetic Hamiltonian that has the usual couplings of 
the ordered double perovskite, and an additional antifer- 
romagnetic coupling when two magnetic B ions are near- 
est neighbour (NN) . The Hamiltonian for the microscopic 
model is: 

H = H loc {rj} + H kin {rj} + H mag {rj}, (1) 

Hioc{ri) = e B Ei VifLfia + e B ' E^ 1 - ViW ia m ia is the 
onsite term where €3 and e#/ are level energies, re- 
spectively, at the B and B' sites. Here / is the elec- 
tron operator referring to the magnetic B site and m 
is that of the non-magnetic B' site. The NN hopping 
term is given by H kin {r]} = -tiJ2<ij>aViVjfLfj^ ~ 
^ £< ii><r (l - ^X 1 - VjW ia m ja - h Y, <ii>a {m + rjj - 

2 r 1i r 1j)(fia m j<r + h' c ')- For simplicity we set all the NN 
hopping amplitudes to be same tx—ti—H = t. The mag- 
netic interaction term consists of the Hund's coupling 
J on B sites, and AFM superexchange coupling J AF 
between two NN magnetic B sites. Thus, H mag {rj} = 

JHiViSi'fL°<xpfip + J AF Y,{ij) ViVjSi-Sj- Here S, is 
the classical core spin on the B site at with \Si\ = 1. 
We take J/t > 1 with J > and J AF \S\ 2 /t = 0.08, 
based on the T^ scale in SrFeOs. We have ignored or- 
bital degeneracy, coulomb effects, e£c, to focus on the es- 
sential magnetic model on the disordered structure. We 
will use a two dimensional model because it already cap- 
tures the qualitative physics while allowing ease of visu- 
alisation and access large system size. The formulation 
readily carries over to three dimensions as well. 
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FIG. 1: (Colour online) First column contains the structural 
motif for four disordered families with progressively increas- 
ing disorder (from top to bottom). We plot the g(ri) = 
(rji — |)e 17r ( Xi+yi ) . We denote the configurations as CI, 
C2, C3, C4 and the corresponding structural order parameter 
has values S = 0.98, 0.88, 0.59, 0.17 from top to bottom. 
Second column shows the ground state spin overlap factor, 
hi = So- Si, where So is the left- lower- corner spin in the lat- 
tice. In the third column, we have shown the corresponding 
NN bond configurations. Here we have three different type 
of NN bonding present between B-B, B'-B' and B-B', repre- 
sented by colours red, blue and green respectively in the plot. 
Lattice size is 40 x 40. 



We have used a real space exact diagonalisation based 
Monte Carlo method involving a traveling cluster approx- 
imation (TCA) 10 to anneal the spin-fermion system to- 
wards its ground state in the disordered background. 

Annealing the electron-spin system down to low tem- 
perature on a given structural motif leads to the mag- 
netic ground states shown in the middle column of Fig.l. 
We plot the spin overlap factor, hi = S .Si, where S 
is the left-lower-corner spin in the lattice. The compari- 
son of the first and second columns in Fig.l indicate that 
the structural and magnetic domains coincide with each 
other. The third column of Fig.l shows the NN struc- 
tural partners. We have three possibilities: B-B, B'-B' 
and B-B', represented by colours red, blue and green re- 
spectively in the plot. 
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FIG. 2: (Colour online) Comparison between the evolution of 
the spin structure factor S(k) at k = (0, 0) with temperature 
for the spin configurations of various disorder families (from 
top to bottom) CI, C2, C3 and C4 obtained from the full 
electronic Hamiltonian with J AF S 2 /t = 0.08 and the effective 
Heisenberg model with Jf/£ = —0.04 and Jaf/£ — 0.065. 
Lattice size is 40 x 40. 



C. Effective Heisenberg Hamiltonian 

Considering the difficulty in doing a spin- wave analysis 
on the full electronic-magnetic Hamiltonian (Eq. [I]), we 
assume that the spin dynamics can be described by an 
effective Heisenberg model 



eff 



{ij} 



(2) 



where {} represents the set of NN and next nearest 
neighbour (NNN) sites. Jij is the effective coupling 
(FM/AFM) between the local moments at and rj sites. 
In our two dimensional ASD configurations Jp operates 
between two local moments when they are at the NNN 
position and J af is active when the moments are at the 
NN position (a B-O-B arrangement). We have estimated 
the effective coupling Jp and J af as follows. For get- 
ting the FM coupling (Jf) we have considered the or- 
dered double perovskite structure. We calculated the or- 
der parameter, i.e, the magnetic structure factor S(k) 
at k = (0,0), as a function of temperature for the full 
electronic Hamiltonian (Eq. [I]) using Monte Carlo sim- 
ulation. We then repeated the same procedure for the 
NNN FM Heisenberg Hamiltonian, defined on only the 
magnetic sites of the double perovskite. We found that 
for J F /t = —0.04, two results matches very well. 

In order to get the AFM coupling we considered the 
ordered perovskite where both the B and B' site carry 
a magnetic moment (mimicking SrFeOs) and computed 
its AFM structure factor peak k=(7r,7r). This model in- 
volves both electronic kinetic energy and Fe-Fe superex- 
change. We find that the result can be modelled via a 
Heisenberg model with Jaf/1 = 0.065. 

Using the couplings inferred from these limiting cases, 
Jp/t = —0.04 and Jaf/1 = 0.065, we studied the bond 
disordered Heisenberg model for the antisite disordered 
DP magnet. We compared the FM structure factor peak 



S(k) at k=(0,0) obtained from the disordered Heisenberg 
model with that from the full electronic Hamiltonian (Eq. 
[T]). The Heisenberg result for the FM structure factor 
S(0,0) as a function of temperature matches very well, 
Fig. 2, with the electronic Hamiltonian result for all ASD 
configurations. This gives us confidence in the usefulness 
of the Heisenberg model for spin dynamics. 



III. SPIN DYNAMICS 
A. Spin- Wave Excitation 

In this section we use the spin rotation technique^ to 
evaluate the spin- wave modes and dynamic structure fac- 
tor at zero temperature. The effective Heisenberg model 
(Eq. |2| can be cast in a form useful for spin wave analy- 
sis by defining a local frame at each site so that the spins 
point along the +z direction in the ground state. We can 
use Si = UiSi, where points along its local z— axis in 
the classical limit. The unitary rotation matrix Ui for 
site Ti is given by 

cos(^)cos(^) cos(Oi) sin(^) — sin(^) 
Ui = -sin(^) cos(^) , (3) 

sin(#i) cos(V^) sin(#i) sin(^) cos(^) 

where and ip are the Euler rotation angles. Now one 
can write the generalized Hamiltonian 

H eff ^2 J ijSi- F ijSji ( 4 ) 

m 

where — UiUj is the overall rotation from one ref- 
erence frame to another and its elements Fff can be 
. . l j 

obtained from Eq. (|3|. 

Applying the approximate Holstein-Primakoff (HP) 
transformation in the large S limit the spin operators 
in the local reference frame become: St - 



2S b u 

S~ = V2S b\ and Sf = S - b\b h where h and b\ are 
the boson (magnon) annihilation and creation operators 
respectively. Only retaining the quadratic terms in b and 
b\ which describe the dynamics of the non-interacting 
magnons and neglect magnon interaction terms of order 
1/5, the generalized Hamiltonian (Eq. [4| reduces to 

« = E^( G <i 6 i 6 J + ^i & *^ +ft - c 0+^( 6 i 6 <+ & 5 6 i)]. (5) 

{ij} 

where Jij = SJij/2, = —SJijFf? and the rotation 

coefficients G? = (Ff?± Fjf) -i( F??^F??). The Hamil- 
tonian (pi) is diagonalized by the transfermation 

bi = V«c n + <*4J, 



(6) 



where S and c are the quasiparticle operators, u and 
v, which satisfy J2 n ( u n u n ~ v n yJ n) = $ij ensuring the 
bosonic character of the quasiparticles are obtained from 



B* 3 



Bi 



4* 



4 



where A %3 = J %3 [G\- + G% ) + eA 3 , ^ = ^ (G?- + <%) 
and = Yljihj + /i*)- Now the spin-spin correlation 
function can be evaluated using the magnon energies and 
wavefunctions obtained from Eq. ([7]), where the excita- 
tion eigenvalues uo n > 0. 



B. Dynamical Structure Factor 

A neutron scattering experiment measures the spin- 
spin correlation function in Fourier and frequency space 
S(k,uj) to describe the spin dynamics of the magnetic 
systems on an atomic scale. From Si = U^Si one can 

express — ^2 U^ a , where a and \i represents the 
y, and z components. Now applying the approximate 
HP transformation to the rotated spins one can write 



Sf =P?& i + gf&t + rf(S-&&), 



(8) 



where /? = +, — and z, and p, g and r are the rotation 
coefficients (given in the Appendix). 

Putting Eq. (J6| in (|8| the space time spin-spin corre- 
lation function can be written as 

S?(t)Sf(0) = 5>±4(*K(0) + B^c m (i)4(0)], (9) 



where the coefficients A and B are expressed in the Ap- 
pendix. In Fourier and frequency space 
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and the total spin-spin correlation function 
S(k,w) = i[5+'-(k,w) + 5-'+(k,w)] + 5 z ' z (k,a;) 



where the coefficient of the delta function 



(ii) 



is the SW weight with B\, = \{B+~ + 5,7+) + Bf? . W£ 

ij ij ij 

is observed as the intensity of magnon spectrum in the 
neutron scattering experiment. 



IV. RESULTS AND DISCUSSION 

We start by presenting the results for magnons in the 
configurations C1-C4 shown in Fig.l and then move to 
an analysis of the linewidth, the estimation of domain 
size, and the contrast with uncorrelated disorder. 



A. Results for AF coupled domains 

Fig. 3 shows the magnon spectra of C1-C4 with ob- 
tained from the Heisenberg model with the FM and AFM 
couplings discussed earlier. In a model with only FM cou- 
plings, i.e., no disorder, we would have obtained only the 
red curve, cj£, for propagating magnons. The striking 
feature in all these panels is how closely the mean energy 
of the magnons follow cj£ despite the large degree of mis- 
location in C2 and C3 and maximal disorder {x ~ 0.5) 
in C4 (refer to the spatial plots in Fig.l). The broad- 
ening, although noticeable in C4, does not obscure the 
basic dispersion. 

Fig. 4 quantifies the mean energy and broadening by 
computing: 



J S(k,uj)ujduj 



5(k, uj)uj 2 duj\ 



■4 



We have shown these two quantities for the C2-C4 struc- 
tures in Fig.l. The Cu\^ have been vertically shifted for 
clarity and the Au^ are superposed as 'error bars' on 
these. It is clear that even in the most disordered sample 
(C4), where the mislocation x ~ 0.4, the broadening is 
only a small fraction of the magnon energy. This will be 
an indicator when we discuss spin waves in an uncorre- 
lated disorder background. 




FIG. 3: (Colour online) Spin- wave spectrum along main sym- 
metry directions of the Brillouin zone for spin configurations 
CI, C2, C3 and C4 (x = 0.01, 0.11, 0.21 and 0.41 respec- 
tively), shown in Fig.l With increasing ASD from CI to C4 
the spectrum becomes broader for a fixed value of momentum 
k. Here J F = -0.04, J A f = 0.065, and lattice size is 40 x 40. 
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FIG. 4: (Colour online) Mean spin wave energy cDk (dots) 
and spin- wave width Ac^k (bars), denned in the text, for the 
correlated antisite configurations C2-C4. The curves are ver- 
tically shifted for clarity. 
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FIG. 6: (Colour online) Mean spin wave energy cDk (dots) 
and spin-wave width Ac^k (bars) for C2-C4 now with Jaf = 0, 
ie, decoupled domains. The curves are vertically shifted for 
clarity. 



B. Broadening: impact of domain size 

There are two ingredients responsible for the spectrum 
that one observes in Fig. 3, (i) the domain structure, and 
(ii) the AF coupling across the domains. To deconvolve 
these effects and have a strategy for inferring domain size 
from neutron data, we studied a situation where we set 
Jaf = in the Heisenberg model denned on the struc- 
tures C1-C4. In that case we will have decoupled FM do- 
mains without any antiparallel spin orientation between 




(JI, -7t) 
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FIG. 5: (Colour online) Spin- wave spectra along main sym- 
metry directions of the Brillouin zone for spin configurations 
CI, C2, C3 and C4 shown in Fig.l (x = 0.01, 0.11, 0.21 and 
0.41 respectively). Increasing fractional weakly coupled do- 
main boundary spins from CI to C4 enhances the SW soften- 
ing near the ZB along [ir, 0] and the spectrum also becomes 
broader for a given k. Here Jf = —0.04, Jaf = 0, and lattice 
size is 40 x 40. 



them. We think this is a interesting scheme to explore 
since the AF bonds are limited to the domain boundaries 
and is not equal to the number of mislocated sites. 

Fig. 5 shows the overall magnon spectrum for this case, 
using the same convention as in Fig. 3, while Fig. 6 quan- 
tifies the mean energy and broadening in this 'decoupled 
domain' case. The absence of Jaf does not seem to make 
a significant difference to the spectrum as a comparison of 
Fig. 4 and Fig. 6 reveal. This correspondence, valid even 
in C4, suggests the following: (i) most of the spectral 
features arise from the domain structure, and the asso- 
ciated confinement of spin waves, rather than the AF 
coupling, and (ii) we can proceed with a much simpler 
modelling of the spectrum and estimation of domain size 
without invoking the complicated BdG formulation that 
AFM coupling requires. 
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FIG. 7: (Colour online) Modelling C2 in terms of a domain 
of size 10 x 10 (left) and of seven domains of size 4x4 (right). 
The corresponding mean energy and broadening are shown 
below. 
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Essentially, much can be learnt from 'tight binding' 
models defined on appropriate stuctures, as happens for 
FM states, without having to invoke the 'pairing' terms 
that arise for AF coupling. A modelling of the full disper- 
sion will require the AF terms as well, but the inference 
about presence of domains, and an estimate of their size, 
need not. We proceed with this next. 

To make some progress in estimating the typical do- 
main size we need a few assumptions; (i) the total degree 
of mislocation, x, should be known, based on the bulk 
magnetisation measurement, (ii) If the overall system 
size in L x L (or equivalent in a 3D model), the number 
of mislocated sites would be xL 2 . (iii) If the domain size 
is Ld then the number of domains within the LxL area is 
Nd ~ xL 2 jl? d . In reality domains need not have one sin- 
gle size, as C2-C4 indicate, but we need the assumption 
to make some headway, (iv) We need to locate these Nd 
domains randomly, in a non overlapping manner, within 
the LxL system, and average the spectrum obtained 
over different realisations of domain location. 

This scheme, carried out for various L^, can be com- 
pared to the full £(k,u;) data to get a feel for the ap- 
propriate Ld- We show the result below for such a tight 
binding exploration for the C2 configuration, modelled in 
terms of different domain distributions that respect the 
same overall mislocation. 

When we compare the ratio of mean broadening to 
bandwidth obtained at different values of Ld (and so Nd) 
with that for the real data, Fig. 4, it turns out that Ld = 
10 provides a best estimate. It also reasonably describes 
the broadening at stronger disorder, C3 and C4, where 
of course Nd is larger. An analytic feel for these results 
can be obtained by considering the modes of a square 
size Ld x Ld under open boundary conditions. 



ground states are disordered ferromagnets but without 
any domain pattern. We computed the magnon line- 
shape in these configurations, and, for illustration, show 
the results for C' 2 and C 3 side by side with their corre- 
lated counterparts C2 and C3. 

There is a striking increase in the magnon line width 
(or Aco>k) in the uncorrelated case. There is almost nine 
fold increase in the magnon line width in C2 and six fold 
in C3 of the uncorrelated disorder with respect to the 
correlated disorder case. 



V. CONCLUSION 

We have studied the dynamical magnetic structure fac- 
tor of a double perovskite system taking into account 
the basic ferromagnetic ordering tendency and the de- 
fect induced local antiferromagnetic correlations. We 
used structural motifs that correspond to correlated dis- 
order, obtained from an annealing process. The results 
on magnon energy and broadening reveal that even at 
very large disorder, the existence of domain like struc- 
ture ensures that the response has a strong similarity to 
the clean case. We tried out a scheme for inferring the 
domain size from the spin wave damping, so that exper- 
imenters can make an estimate of domains without hav- 
ing spatial data, and find it to be reasonably successful. 
We also highlight how the common assumption about 
random antisites, that is widely used in modelling these 
materials, would lead to a gross overestimate of magnon 
damping. In summary, dynamical neutron scattering can 
be a direct probe of the unusual ferromagnetic state in 
these materials and confirm the presence of correlated 
antisites. 



C. Contrast with uncorrelated antisites 

In modelling the antisite disorder much of the earlier 
work in the field assume the defect locations to be ran- 
dom. We have followed the experimentally motivated 
path which suggests that the mislocated sites themselves 
form an ordered structure separated from the parent (or 
majority) by an antiphase boundary. The sources of scat- 
tering are the boundary between these domains rather 
than random point defects. Since much of double per- 
ovskite modelling has assumed the random antisite situ- 
ation it is worth exploring the differences in the magnon 
spectrum between correlated and uncorrelated antisites. 

We have already seen the results for correlated disor- 
der for different degrees of mislocation, x. We generated 
uncorrelated antisite configurations with the same x by 
starting with ordered configurations and randomly ex- 
changing B and B' till the desired degree of disorder is 
reached. These configurations naturally do not have any 
structural domains. Annealing the full electronic Hamil- 
tonian on these configurations, call them C(, .., etc, 
down to low T, leads to the magnetic ground states. The 
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VII. APPENDIX 

The rotation coefficients are 

pt = yf§{Ur±U? v )-i{U? x TUr v ) 

rf = Uf x ± iU* y (12) 
pi = Uf z - iUf 

q? = ur+iu? z 

r z = TJ ZZ 

1 i ^ i 
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FIG. 8: The left set of panels correspond to mislocation x = 0.11 where we compare the magnon spectrum for uncorrelated 
disorder (left) with correlated disorder (C2) right, the top panels refer to the structural pattern, the middle to the magnetic 
ground state, and the bottom to the magnon response. The right set of panels refer to x = 0.21, and the same indicators as for 
the left panels. Notice the remarkably broader lineshape for the uncorrelated disorder case where it is difficult to make much 
of a correspondence with the clean dispersion. 



And the structure factor coefficients are 

ij ±i 1 j 1 j 1 ± j 1 j 1 ± j 1 j 
+q?q^u 7 i n *v™ - S x rfr^(u^uf + uf* ' v%) 
B± = q?pPv?vf +pfqfuY l uf +pTp$u?vf 
+q?q^v™uf - S x rfr^v^vf + vfvf). 
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